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corrections due to instanton sectors k = dbl, contributing to the propagator, 
Q ■ are shown to be just the homogenous terms admitted by the Dyson-Schwinger 

^ | ■ equation for S. In the case of the 4-fermion function also sectors k = ±2 are 

t^j- , included into consideration. The quark condensates are then calculated and 

(N ■ are shown to satisfy cluster property. The ^-dependence exhibited by the 

I Green functions corresponds to and may be removed by performing certain 

!>■ ' chiral gauge transformation. 

O ! 11.10.-z;11.30.Rd;12.90.+b 
(N . 



00 
On 




> 

X 



Typeset using REVTeX 



*Elecronic mail: torado@fuw.edu.pl 



1 



I. INTRODUCTION 



Many phenomena in present-day theory of elementary particles can be only described 
in the nonperturbative language. In the first place (though not exclusively) it refers to 
the theory of strong interactions — quantum chromodynamics — where the only existence 
of mesons and baryons requires allowing for such effects as quark and gluon confinement. 
However QCD, based on colour SU(3) symmetry group, introduces a variety of coupled 
fields, which makes the perturbative calculation very complicated, and nonperturbative one 
extremely perplexed. 

Fortunately we have at our disposal a couple of examples of model theories, that, thanks 
to their mathematical simplicity, allow for investigating certain nontrivial and nonperturba- 
tive aspects of more complicated phenomena. Among most fertile ones one should mention 
the so called Schwinger Model — the electrodynamics of massless fermions that, in this work, 
will be called 'quarks', in two space-time dimensions. This model was originally proposed |TJ 
as an example of a field theory in which the existence of the local gauge symmetry does not 
entail the zero mass gauge boson — a fact which should not be pushed aside in electroweak 
interactions theory. It turned out soon, however, that it has many other nontrivial proper- 
ties, which can be interesting from the point of view of both, strong and weak interactions. 
Above all one should mention here the confinement of fundamental fermions and the axial 
anomaly. 

The similarity between quantum chromodynamics and the Schwinger Model appears also 
(and this will be particularly essential for this work) in the presence of topological effects: 
the existence of instanton sectors leading to the emerging of vacuum and nonzero 'quark' 
condensate. 

It is a known fact || that if the theory contains massless fermions (at least one) the am- 
plitude of the tunnelling transition between different topological vacua is diminished to zero. 
Mathematically it is expressed through zero value of the Euclidean Dirac operator determi- 
nant: det[i ^-e|], because of the appearance of the zero eigenvalues and eigenvectors when 
A bears the instanton index 0-0]. The functional integral over fermionic degrees of freedom, 
corresponding to the transition amplitude in question, is then zero and the tunnelling phe- 
nomenon disappears. On the other hand, however, the notion of the ^-vacuum does not lose 
its sense since topological vacua do not exhibit the so called cluster property M. The reason 
is that, despite the appearance of the anomaly, a conserved chiral charge still can be defined, 
the topological vacua being its eigenstates with different eigenvalues. The matrix elements 
< m\Q(x)\n > of certain local operator Q(x) between such vacua (|n > and \m >) are then 
nonvanishing if the operator Q changes the chirality (the Hamiltonian does not belong to 
that type of operators and hence the above mentioned lack of tunnelling). If one of the 
topological vacua were taken as the true vacuum, for such a type of operators Q one would 
obtain automatically < Q(x) > vac = 0. A product of such operators P(x,y) = Qi(x)Q2(y) 
can, however, have the nonzero vacuum expectation value, for instance if it is constructed 
in a chiral invariant way (e.g. if Qi changes the chirality by k and Q2 by —k). Now the 
requirement of clusterization in the limit \x — y\ — > 00 leads to contradiction || 

^< P > 

vac < Qi > vac < Q2 >vac= . 

If we then have to do with vacuum expectation values for chirality nonconserving operators, 
we have to include in the calculation different topological vacua. It is just that category that 
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fermion, and fermion-boson Green functions belong to. Because of the 2 : 1 correspondence 
between the chiral charge and the topological index of a vacuum, the objects bilinear in 
fields \I/ (propagators) require the inclusion of the instanton sectors 0, ±1, and 4-fermion 
functions — sectors 0, ±1, ±2. 

All these facts are, of course, well known and applied successfully to the calculations of 
the quark condensate < ^(x)^(x) > vac in the Schwinger Model, as well as to the verification 
of the revival of the cluster property for < \I> (x)^ (y)^ (y) > vac 

The topological methods, however, have not been till now applied to the calculation 
of the full 2- and 4-point Green functions and not only condensates. This computation 
will constitute the main goal in the present work. The paper is organized in the following 
way. In Section |1| we concentrate on the quark propagator. From the considerations of the 
Dyson-Schwinger equations one knows [p~3| , [T^[] that, beside the famous Schwinger solution U 
satisfying {S, 7 5 } = i.e. S ~ 7 /i (...) At , they admit also other terms, let us call them S', 
for which [S",7 5 ] = 0. That means that S' ~ 7 5 (---) + The authors raised then the 



question concerning the interpretation of such terms. Extending slightly the analysis of [14 



we show, in Section |11 A| , which the most general form of the quark propagator, resulting 
from the Dyson-Schwinger equations, is. Next, in Section |II B| using the methods of |10 



applied there to find the value of condensate, we calculate the full propagator and show 
that the additional terms 5" are just those emerging from nonzero instanton sectors. In 
section |T| we perform the similar calculation for the 4-fermion function. The analysis is 
here more complicated since it now involves the 0, ±1, ±2 instanton sectors and also because 
of the tensor structure of G, which bears now the four spinor indices. In the final section we 
verify the cluster property and have a look on the 9 dependence of the full Green functions. 
As one knows, in massless theory the parameter 9 may be removed by the convenient chiral 
gauge transformation. This conclusion will find its full confirmation in our expressions for 
the Green functions. 



II. INSTANTON CONTRIBUTIONS TO THE QUARK PROPAGATOR 

In this section we would like to concentrate on the fermion propagator. Being defined as 
the vacuum expectation value of the product of \1/ and \I> it should acquire additional terms, 
beside that found already by Schwinger, originating from instanton sectors ±1. 

We start with summarizing briefly the conventions used in this work. The two- 
dimensional Lagrangian density of the Schwinger Model with the gauge fixing term has 
the following form 

C(x) = (z 7 ^ - erfA^x)) ¥(x) - h^(x)F flu (x) - ± (d.A^x)) 2 , (1) 
where for Dirac matrices 7 M we assume the representation in which all gamma's are real 

Mr;). '-Mi-!)- (2) 

The metric tensor g^ v and the totally antisymmetric symbol e^ u are defined as follows 
g ™ = -g* = 1 , £ ^ = -s W = 1 , e 00 = e 11 = . 
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A. General considerations 



Having specified the Lagrangian we can derive, in the standard way, the Dyson-Schwinger 
equations for the propagators for the two basic fields in the theory. It is well known 



that the result for the gluon Green function is simply that it acquires a mass equal to y e 2 /V- 
We recall that in two dimensions e is a dimensionful constant. No other contributions from 
the nonzero instanton sectors come into play since A 11 does not change the chiral charge, 
similarly as interactions with photons do not change the electric charge (which is no more 
true in nonabelian theories). We, therefore, take the full, dressed gluon propagator as already 
known to be (in momentum space) 



D^{k) 
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1 k^k 1 



k 2 J e 2 /vr - k 2 X (k 2 ) 



2\2 ' 



(3) 



and concentrate only on the quark sector. The Dyson-Schwinger equation, as always, chains 
up the 2- and 3-point functions which recurrence (continued to infinity) usually efficiently 
precludes us from solving it 



jS(p) = 1 + zeV / ^S(p + k)r v (p + k,p)S(p)D^(k) . 



(2k) 



(4) 



The prominent and well known advantage of the Schwinger Model is that its Lagrangian is 
invariant under two types of gauge transformations: ordinary and chiral ones. This leads to 
two kinds of Ward identities [ |H]| , p] -|TP|] which relate the projections of r M (p + fc,p) on k^ and 
on e^ v k v with quark propagator. In two-dimensional world they are sufficient to reconstruct 
the full vertex Y^ 1 and decouple the equation from the infinite hierarchy 



s{p + fc)r*(p + k, k)s( P ) = K - [s(p) - s( P + k)\ - 



-y 5 s(p) + s(p + k)r 



(5) 



If we now adopt the Landau gauge (A — > oo) we see that the longitudinal part of @ does 
not contribute since D^ u {k) becomes purely transverse and equation (|4]) may be given the 
following closed form 



ffS(jp) 



le 



d 2 k 



(2w) 2 k 2 (e 2 /K - k 2 -ie) 



i 5 s{p) + s{ P + k)r 



(6) 



where the relation e Ma 7 M /c a = /^7 5 has been used. In the co-ordinate space there is a known 
factorization and the equation is simplified to 



i @ x S{x) = 5^ 2 \x) - e 2 [fcP(x)] i 5 S(xh 5 
The function (3 is here defined as follows 

1 




(jp. _ e 2 _|_ 2e)(p 2 -f ie) 



(7) 



(8) 



-f + Ie + In Je 2 x 2 /AK + fH^Ue 2 x 2 /ir 



7s + In J—e 2 x 2 /4ir + K (J — e 2 x 2 /k 



x timelike 
x spacelike 
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and is, in fact, a function of x 2 only. Symbol 7^ denotes the Euler constant and functions 
Hq 1 ^ and K are Hankel function of the first kind, and Basset function respectively [20 



Assume now the most general matrical structure that is admitted in two dimensions, in 
co-ordinate space, for the fermion propagator S 

S(x) = S (x)A(x 2 ) + B(x 2 ) + 7 5 C(x 2 ) + 7 5 tfD(x 2 ) . (9) 

The term proportional to ft has been chosen, for convenience, as explicitly containing the 
free propagator Sq(x) = — ^ x /_ i£ - It is, of course, only a question of redefining the co- 
efficient function A(x 2 ). In what follows we will omit epsilons specifying that we have to 
do with causal propagator. We now insert (Q) into the co-ordinate space Dyson-Schwinger 
equation ([7p and exploit the fact that for all unknown functions, as well as for (3 function, 
we can write j$ x F(x 2 ) = 2 ftdF/dx 2 = 2 $F'(x 2 ). This allows us to deduce the set of four 
differential equations for four functions to be found. This set arises if one takes the trace 
of (0) with successive insertion of 11, 7 5 , and 7 s on both sides. These are simple first 
order equations 



A' = -ie 2 (3'A , B' = ie 2 f3'B , C = ie 2 (3'C , D' 



\ + ie 2 (3' 



x 2 



D , (10) 



with the initial condition A(0) = 1 originating from the cancellation of the Dirac delta 
functions in (0). This set of equations may immediately be solved and we obtain the most 
general form of the propagator that is accepted by the Dyson-Schwinger equation 



S{x) = S {x)e' ie2 ^ x) + B e ie2 ^ x) + C ^ 5 e ie2 ^ x) + B^-^-e^^ , (11) 



where constants B , Co, -Do remain unknown. 



In what was stated above we did not move to far from what had been done in [14|, except 



that we found two additional terms [C and D) in the most general structure of S. In what 



follows, however, we will show that the terms with unknown constants in ([11]), except for 
the last one, arise as a result of instanton effects. The first term is just the well known 
Schwinger solution. 



B. Explicit calculation of the quark propagator 

The instanton contributions to the quark condensate were calculated in a systematic way 
in ||1Q|| . Below we extend this calculation and find the expression for the full propagator. 
We begin with substituting into the generating functional, defined as usually as 

Z[ v ,rj, J}= J vqvqVAe'f^l^^^M , (12) 

the following form of the gauge potential 

A*{x) = A {0) »(x) + e^d v b{x) , (13) 

where A^ ^ is certain new potential satisfying some special conditions at space-time infinity 
(this point will be discussed afterwards) and b is the external scalar function to be chosen 
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later for our convenience. Since ([H|) constitutes a simple shift we can now easily pass from 
the functional integral over A to that over A^°\ It is known that the coupling term e^> A^'ty 
may be gauged away if we introduce new fermion fields defined by the relations 



= W (x)e ier ^ x S d 2 zA(x-z)A { °\z) dcf ij>{x,A {0) ) ? 



(14) 
(15) 



with A(x — z) satisfying: D x A (x — z) = 5^ 2 \x — z). The above gauge transformation is 
an element of U(l) <g> {/a(1) group. While the Lagrangian ([!]) is invariant (apart from the 
gauge fixing term) under that kind of transformations, the fermion path integral measure 
in (12) is not |2l|-p3|. This anomalous behaviour generates a mass term for the gauge boson. 
After this transformation is performed in (|I2|), the dependence reappears in the source 
terms through the functions 4>(x,A^) and (f)(x, A^ ') defined by equations (|14]JT5|) . The 
substitution A^ — > ^ allows us to perform the remaining gaussian integral over A^ [24 
and we obtain the following expression for Z 



Z[t), rj, J] = N b exp < % \ d 2 x 



-bn 2 b + e^Jfivb 



exp < % 



d x 



ma 



ee 



x 



(16) 



exp |-- J d 2 xd 2 y [j» + (□ + e 2 /n)e fia d a b\ A Mi/ (x - y, e 2 /n) [j v + (□ + e 2 /n)e uf3 d p b 

where Nf, is certain constant (b stands for boson). The massive propagator A^(x — y, e 2 /ii) 
satisfies the equation 



□ + 




d d \ 



□ 



pS-'V 



A- ( x -y,e 2 /7r) = 5^(x-y)gl 



(17) 



The main point is now the evaluation of the fermion path integral in (|T^) . Let us denote 
by X the following expression 



X = V^Vm exp \i / d 2 x V (i ft - ee^>y^d u b) * + ffyff + #7/ 



(18) 



where primes are used to avoid writing explicitly the factors e - ^ multiplying the external 
sources. This expression is naturally proportional to the determinant of the Dirac operator. 
From (|18|) it is evident, however, that it strongly depends on the choice of hitherto undefined 
function b. It turns out that for certain choices of b it may even vanish so one could proceed 
here with caution. First of all we temporarily go over to the Euclidean space since the 
properties of the Dirac operator are there mathematically more rigorous. The transition to 
this space is defined by the following substitutions: x° — > — 1x2, A — > iA 2 , d° — > 182- The 
Euclidean metric tensor g^ u is defined as and for e^ u we have E\ 2 = —£21 = — i- In what 
follows the same symbols as before will be used for denoting the Euclidean counterparts of 
earlier defined quantities. We hope that it will not cause any confusion and the passing back 
to the Minkowski space will explicitly be stated. Before we calculate the quantity flTB] ) we 
have to make some remarks on the instanton sectors. 
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Consider the most general form of the field A in the 2-dimensional world 



A^x) = d^a(x) - e^ v d v b{x) , (19) 
The so called Pontryagin index for the gluon field may be defined as 

16- f 16- f 

U =^J d2xe af3F af3 = - — J d 2 xU x b( X ) . (20) 

Due to the Euclidean nature of the space-time the d'Alambert operator n x is here, naturally, 
the same as the Laplace operator. Assume now that we take the function b(x) in certain 
specific form 



h W {x) = ±k1n(*±£) , (21) 



where k is an integer number. After evaluation of (|20| ) we immediately obtain v — k i.e. 
the field A^p defined as £ a p^^ bears the index v — k. Since the Pontryagin index (ED) is 
linear in gluonic field the total index of the sum 

A<g>+e als dpbW , (22) 

is equal to k, and thanks to the complete freedom while choosing A^ (which is restricted to 
the sector 0) constitutes the most general form of the field in the k instanton sector: 
j±(k) represents a path (in the sense of Feynman path integal) connecting two distinct and 
topologically inequivalent vacua for t = ±oo. The true vacuum of the theory, the so called 
9- vacuum, is now the superposition of these topological vacua 

oo 

\e >= eine \ n > > ( 23 ) 

n=— oo 

and the generating functional calculated in this new vacuum has the form 

oo 

Z[ V ,%J]= £ e^zWfo^J], (24) 

k=— oo 

where the summation runs over instanton indices rather than over topological vacua ones. In 
each sector the appropriate is calculated with the restriction on the values of the vector 
potential to those defined by (|2"2"D. Now we are ready for considering the contributions from 
separate terms of the sum in (|24]) — the different instanton sectors. In compliance with 
what was said in the Introduction about chirality conserving and nonconserving operators 
and their matrix elements taken between different topological vacua, we have Z[0, 0, 0] = 
Z(°)[0, 0, 0] = N f ■ N b and for the propagator we can immediately write 



k^oc z{ ) [°. °. °] fc^oo 6r)(x)5r)(y) 



■ (25) 

J,r],rj=0 



Similar formulae are valid also for other Green functions and will be exploited in the following 
section while dealing with 4-point functions. The sum in (p5|) formally extends to +oo but 
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practically depends on the properties of the operator the vacuum expectation value of which 
is being considered. In the case of a propagator this operator is simply a product and 
the whole sum reduces to three terms corresponding to k = 0, k — +1 and k = — 1 which 
means that chiral charge of the contributing vacua can, at most, differ by 2. It is not difficult 
to observe that the case k = corresponds to the well known Schwinger solution since it 
simply requires putting = 0, and is represented by the first term of (|TT|). The nontrivial 
topological effects come from two other sectors. We concentrate below on the sector k — +1, 
the calculation for k — — 1 being analogous. The applied method is here that of |L0| and 
therefore we do not plunge into details and point out only the main steps. 

For k = 1 we have to substitute b = b^ 1 ' into (p~8|) (or rather into its Euclidean version) 



and perform the \1/ and \1/ functional integral. By the Atiyah-Singer index theorem ||-[7j 
the massless Dirac operator in the external field bearing the nonzero Pontryagin number 
possesses zero modes. For k = ±1 there is only such mode and, for gamma matrices 
conventions defined by (0), it has the form (for the discussion and construction of the zero 
modes in the Schwinger Model see |10| , [25| J 

, , 1 



Xo{x) 



2tt 
1 



X 2 + x 2 
1 

27 \X 2 + x 2 



1/2 



1/2 



for k = +1 



for k 



(26) 



(27) 



In the 2-instanton sector, which will be dealt with in Section |TTTj, two zero modes appear 
1 A 3 / 2 f0\ , , 1 . x \ 1/2 f0 



/2^X 2 + x 2 \ l 
1 A 3/2 / 1 
/2^X 2 + x 2 [ 



Xi{x) 



2vr 
1 



Xi - IX 2 



xi + ix 2 




The whole set of eigenmodes will be denoted by \n- We have 



a n \n , * = b nXr. 



for k = +2 , (28) 
for k — —2 , (29) 

(30) 



where a n and b n are the Grassmann coefficients. We now pass from the integration over \l/ 
and \l/ to that over a n and b n . The presence of a zero mode means that neither do nor bo 
appear in the bilinear term in the exponent of (|18D and only source terms in the integration 
over these two parameters are involved. The rules for the Grassmann integration demand, 
however, that from the whole exponent only the linear part contribute. The integration over 



nonzero modes may be done in usual way, the details of which can be found in ||10|| . In that 
way we arrive at the following expression for the quantity X 



X = N f \[ d 2 xd 2 yf]'(x)xo(x)xUyW(y) exp - / d 2 xd 2 yrf{x)e- i ^- £ ^ bW) 



S (x-y)e i ^ y '- £ ^ bW) r ] , (y) 



exp 



_ d 2 x—bWnbW 

2tt 



(31) 



So being the Euclidean free fermion propagator: i x S o (x — y) = —5^(x—y). The additional 
factor A arises from the change of variables ( |30"D since the dimension of field \l/ in the Feynman 
path integral measure is 1/2 and that of ao and bo — 0. 



S 



Now, when the integration over fermionic degrees of freedom has been performed, and 
the properties of the Euclidean Dirac operator have been exploited, it is convenient to come 
back to Minkowski space, where the final formula for the propagator, analogous to flTT|), is to 
be obtained. It should, however, be emphasized that whenever required by the mathematical 
rigour, the corresponding expressions in Euclidean space are presumably considered. 

The differentiation over external currents 7] and rj, as required by (^), leads to a very 
simple expression if one keeps in mind that finally we have to set all the external sources to 
zero 



6 2 X 



5r](x)5rj(y) 



. iV/ A e -^(^' 5 / i ' 5J ) Xo (a;)x^(2/)e i ^^ fi / iW ) e i / d2a: ^ 6(1)D6Cl) 



(32) 



77,77=0 



From the definition (|26| ) of the zero mode it is evident that the product X0X0 has the matrical 
form I (1 — 7 5 ). The functional (ft has similar structure (a linear combination of 11 and 7 s ) 
and therefore both quantities commute with each other which leads to the considering of 
the operator 



exp 



-i(ft(x, S/iSJ) + i(ft(y, 5/iSJ) 



(33) 



where functional differentiations act on the J dependent argument of the last exponent 
function in ([16]). Let us now introduce the (nonconserved) current AT M by the relation 



{x,A)-ij>{y,A) 



d 2 zA fJl (z)K, IM (z;x, y) 



and satisfying 



W x,y) = e [S {2) (x -z)- 5^{y - z) 



(34) 



(35) 



which is evident by virtue of the definitions ( ]T^JT5| ) of (ft and (ft. In this notation the op- 
erator (P^) is just: exp [/ d 2 zJC^(z; x,y)5/5J^(z)} and simply shifts the argument of the J 
dependent functional 



exp 



J d 2 zW{z;-, 
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5J»(z) 



F[J}=F[J + IC} 



(36) 



Having found expression ([32]) for the second derivative of X, and collecting (|16| ) together 
with fl2"5]) we are in a position to write down the formula for the one-instanton contribution 
to the quark propagator 



S^(x -y)= i\e l9 Xo (x)xUy) ' exp {- / d 2 u [&«(«)□ (□ + e 2 /n) b^(u) 



(37) 



d ud w 



/C M (w; x, y) + (□ + e 2 /vr) d a b {1) (u)l A M „ {u - w) \K v (w\ x, y) + 



^ (n + e 2 /n)dpb {1 \w 



In the exponent of ( p7|) all functions are perfectly known, so it is not difficult, although 
lengthy, to evaluate all expressions. We skip this calculation here and only note that all 
terms may be divided onto three groups: those not containing current /C M , linear, and 
finally quadratic in /C M . 
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1. It is not difficult to check that both terms which do not depend on /C cancel each 
other. 

2. The terms linear in K, are strongly simplified if one notices that thanks to the pres- 
ence of e^ u from the whole expression for only that proportional to metric tensor 
contributes. If one now exploit the known identity for gamma matrices in two dimen- 
sions 7^7 5 = e^Tn and observe that finally 7 s may be replaced with -1, because of 
the structure of the matrical coefficient (xoXo ) i n (pTD, we see that the total effect of 
these terms is just to cancel X0X0 ( U P to together with the coefficient A. This was, 
naturally, expected since the conclusive formula should not contain any trace of the 
particular choice of 

3. The calculation of the term quadratic in K, is also elementary and the result turns out 
to be 

lE + \\n-^ + ie 2 (3{x-y) . (38) 



This allows us to complete the formula ( j3"7| ) for S^. The twin calculation in the k = — 1 
instanton sector gives the analogous outcome with the reversed signs of parameter 9 and 
matrix 7 5 . Finally, the quark propagator with contributions from all sectors is 

S{x) = S (x)e- ie2 P {x) + (cos 6 - 2 7 5 sin#) e 7B+ie2/3(x) . (39) 

This proves that the second and the third term in (0), admitted by the Dyson-Schwinger 
equation, do actually materialize as a result of instanton effects. The last term does not 
reappear in (|39"D, and in fact may not do, by virtue of simple dimensional consideration. 
From (|TTD it is obvious that the constant D has to be dimensionless and simultaneously 
it must tend to zero when e 2 — > (the free case). Naturally no such constant exists in a 
world in which the only dimensionful constant is just e 2 . From similar dimensional analy- 
sis it is obvious that constants B and Cq in ( ]TT| ) must be linear in e which is in fact the 
case in (|39|) . One should note in this point that instanton effects usully manifest them- 
selves in a nonperturbative way, for instance through the appearance of negative powers of 
the coupling constant. As an example we can quote the one-instanton tunelling amplitude 
in four-dimensional YM: < 0|1 >~ exp(— 8n 2 /g 2 ) [2~7|. In the present case, however, we 



cannot construct e 2 -dependent dimensionless constants, and the intuition concerning non- 
perturbative contributions should be confronted with dimensional considerations. In the 
last section we come back to the formula (|39|) and analyze it from the point of view of the 
quark condensate. 

It should be also observed that other functions with two quark 'legs', as for instance 
the vertex function, are still completely expressible through the propagator S as discussed 



in |L9| and the relation 

[ST^ (x, y; z) = iS^(x - y)Y $ z A(y-z)—i ^ A (x - z)S^ k \x - y) , (40) 
as well as its counterparts for higher functions, hold separately in each instanton sector. 
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III. INSTANTON CONTRIBUTIONS TO THE 4-FERMION FUNCTION 



In compliance with our former discussion, the expectation value of the product 
acquires contributions from up to k = ±2 instanton sectors. As in the case of the propagator, 
the 4-fermion Green function in a certain specific sector is defined by the functional derivative 
of the generating functional 



ike 



8 4 Z^[r],fJ, J] 



Z[0, 0, 0] S7] a (xi)Sr) b (x 2 )5ri c (x 3 )5r] d (x4) 



(41) 



r],rj,J=0 



For the sector k = the appropriate expression for the Green function was found in our 
previous work We, therefore, begin with considering the case k = 1 recalling here only 
the final result for G^ 



G < S, cd (x 1 ,x 2 ;x 3 ,x 4 ) 



<S()ac(xi ~ X 3 )S bd(x 2 ~ X 4 ) + (S Q (xi - ^3)7 5 ) ac (<So(x 2 ~ X A )^ 



bd 



exp 
1 

+ 2 
exp 



ie 2 ((3(xi - x 2 ) - f3(xi - x 3 ) - j3(x 2 - x 4 ) - (3(xi - x&) - (3(x 2 - x 3 ) + (3(x 3 - x 4 )) 



<Soac(xi - X 3 )S 0bd (x 2 - X 4 ) - (S (Xi - X 3 )^f 5 ) (S (x 2 - X A )"f' . 

\ / ac \ / bd 

- ie 2 ([3(xi - x 2 ) - (3(xx - x 3 ) - /3(x 2 - x A ) - f3(xi - x 4 ) - (3(x 2 - x 3 ) + [3(x 3 - x A )) 



-(c, x 3 



d, xA 



(42) 



The evaluation of has already been done in the previous section so we are able to 
immediately write down the following equation for G^ 1 ' 



G ab,cd( x ^ X ^ x 3, x 4 ) = i\e te 



{a,x x <r+ b,x 2 ) 



exp 



J (fxb^D (□ + e 2 /7r) d 2 xd 2 y (>(x) + (□ + e 2 /n) d a b w (c 



A M „(x - y) (r(y) + e** (□ + e 2 /n) dpb^(y) 



(43) 



taken at J = 0. We now recall the definition (|34T) of the current /C from the preceding 
section and J introduced in |T^] by the relation 



<f>(x, A) - </>{y, A) = -J d 2 zA\z)J^z-x,y) . 

The both currents satisfy the equation 

JC^z; x, y) = J^{z; x, y) - 2e 7 V"^ A {x - z) . 
Using this notation we can rewrite ([13|) in the form 



(44) 



(45) 
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G ( * 1) (x 1 ,x 2 ;x 3 ,x 4 ) = i\e l9 (xo(xi)Xo( x 3) ® Sofa - X4)) ^exp J d 2 zlC fl (z;xi,x 3 )j-^ 



) exp 



cxp 



% J d 2 zj ll {z;x 2 ,x 4 )e^d u b {1 \z) + J d 2 zj^{z;x 2 ,x 4 ) — 
J d 2 xb {l) U (□ + e 2 /vr) d 2 xd 2 y (>(x) + (□ + e 2 /ir) d a b {l \: 



X) X 



/\^(x-y)(J u (y)+s u(3 (n + e 2 /it)dpb {1) (y)^ + antisymmetrization , (46) 

where the 'antisymmetrization' is defined by substitutions in (|43"1). The obvious matrical 
notation has been introduced above to simplify the expression and avoid explicitly writing 
four spinor indices. (One should keep in mind that both currents K, and J are matrices in 
spinor space.) The functional derivatives over current J, similarly as in the expressions of 
the previous section, lead to shifts of J and we obtain 

G») (!tl ^;x, t x 4 ) -i^LtehJfe) - *.)«p (< JfzJAr,*,,*.) x (47) 



z^d v b {1 \z] 



■ J d 2 xb {1) U (□ + e 2 /n) - l -J d 2 xd 2 y{p i {x- x u x 3 ) ® 31 
1 ® J(x; ar 2 , x 4 ) + ^ Q (□ + e 2 /rt) d a b {l) (x)) (x - y) (/C"(y; a*, x 3 ) ® U + 
1 <S> J(y, x 2 , x 4 ) + £ w/3 (n + e 2 /V) dpb (1 \y)^j + antisymmetrization , 



+ 



In the last exponent of (f47|) we recognize the expressions similar to those of fl37j) although 
the tensor structure is now much more perplexed. Nevertheless all functions are known and 
the evaluation of ( f4?j) is only a matter of patience. We do not intend to go into details and 
only indicate the main points. 

1. The term quadratic in K is just that of the previous section multiplied (tensor product) 

by a. 



2. The term quadratic in J was calculated in our previous work \W 

3. The two terms containing squares of b^ cancel each other. 

4. The terms containing products of J and cancel with the first exponent in (fi7|). 

5. The terms containing products of K, and b^ cancel with the appropriate part of the 
coefficient Xo{ x )Xo(y) ( see Section |IIB|) . 



6. The product terms of both currents K, and J may be evaluated in the straightforward 
way and we obtain 



ie 7 5 <g> 7 5 [(3(x 1 - x A ) - (3{xi - x 2 ) - f3(x 2 - x 3 ) + (3(x 3 - x 4 )] . 



(48) 



With all above taken into account we can now the arising exponent in ([IT]) give the following 
simple form 
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cxp 



all <g> 11 + &7 5 ® 7 5 



i (n ® 1 - 7 5 ® 7 5 ) exp(a - 6) + 
- (l ® 1 + 7 5 ® 7 5 ) exp(a + 6) , 



(49) 



where a and 6 are certain functions expressible by combinations of /?'s. 

The calculations for k — — 1, which are practically identical, allow us to write down the 
final expression for the one instanton contributions to G 



G (+1) (xi,x 2 ;x3,x 4 ) + G ( 1) (xi,x 2 ;x 3 ,x 4 ) 



, -e 7B (cos 9 — i'y 5 sin 9) ® iS (x 2 — x 4 ) x 

^11 ® 11 — 7 5 <E> 7 5 ^) e ie2 [/ 3 ( x i- x i)-/ 3 ( x 2-X4)-/3{xi-X4)+P(x 1 -X2)+l3(x2-X3)-f3(x 3 -X4)] _|_ ^gg^ 

^11 <S> 1 + 7 5 ® 7 5 ) e ie2 lP( x i- X: i)~P( x 2- x i) + P( x i- x 't)-P( x i- x 2)-P(x2-x 3 )+f3(x 3 -x i 

+ antisymmetrization . 

In the sector k = 2 we use again the formula fl¥ip . There are now two zero modes. 
This means that after functional integration over Grassmann variables, analogous to that 
performed to obtain flST]), the product of two 77's and two 77's appears. These four sources 
have to saturate all four derivatives in (f41|), since otherwise we would get zero after setting 
77,77, J = 0. That in turn means that, in the formula for G^ 2 \ the current J will not appear 
and 



G (2) (xi, x 2 ; x 3 , x 4 ) = -\ 2 e 2ie X o{xi)xt(x3) ® Xi(^)Xi M exp - J d 2 xb {2) U (□ + e 2 /7r) 6^ 
y d 2 xd 2 y(lC f "{x; x x , x 3 ) ® 11 + 1 ® /C(x; x 2 , x 4 ) + e Ma (□ + e 2 /^) d Q 6 (2) ( 

A^(x - y) (/Cfo; X!, x 3 ) ® 1 + 1 ® /C(y; x 2 , x 4 ) + (□ + e 2 /vr) ^6 (2) (y) 
+ antisymmetrization . 



(51) 



We are not going to repeat below the steps leading to the final result, since the modification 
of previous formulae is only slight. The 2-instanton configurations contribute then in the 
following way 



G^ +2 \x 1 ,x 2 ]x 3 ,x 4 ) + G ( 2) (xi, x 2 ; x 3 , x 4 ) 



256vr 4 



M (a - 7 5 ) ® (a - 7 5 ) (*§ + 4) 

^a + 7 5 ) <g> (a + 7 5 ) (-ar§ + :4)(x 4 + x 4 ) exp {ze 2 [/3(xi - ar 4 ) + (52) 
/3(x 2 — x 3 ) + f3(xi — x 2 ) + f3(x 3 — X4) + f3(xi — x 3 ) + (5{x 2 — x 4 )]| + antisymmetrization . 



~x 4 ~j~ x 4 ) ~\~ c 



If we now expand the tensor notation of ( |50| ) and ( |52"D into explicit spinor indices and perform 
the full antisymmetrization (which also restores the apparently broken Lorentz invariance) 
as defined by fl43|) we can collect together all the contributions 



G a b,cd{xi,x 2 ; x 3 , X4) — G^ cd (xi,x 2 ; x 3 , x 4 ) + G a ^J d (xi,x 2 ; x 3 , x 4 ) + G ab J d (xi,x 2 ; x 3 , x 4 ) + 

Gdb%(x u x 2 ; x 3 , x 4 ) + G ( ri%(xi, x 2 ; x 3 , x 4 ) . (53) 
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Since we now dispose the complete formula for the 4-fermion function, we can extend 
the analysis carried out in [lSj], regarding the possible existence of the Schwinger pole in 
the t-channel (the quark- ant iquark scattering), over contributions from higher topological 
sectors. To this goal we identify the appropriate co-ordinates of incoming and outgoing 

X4 and perform the Fourier transform of the expression in 



particles: u = f X\ 



def 

= x 3 , v = x 2 

the new variable z = f v — u. This identification (which should be treated as an appropriate 
limit) has to be done in a symmetrical way as described in the quoted paper. After these 



manipulations we see that the formula (45) from |T9[ acquires additional terms from the 
k 7^ instanton sectors 



G ' polar (k*) 



1\2 



ie 



47T 7 °® 7 V 2 -e 2 /7r) + 47r 3 /2 e7£ 



(7 ® 7 5 e"^ 5 



k 



ie 



(k 2 -e 2 /n) 8tt 



II ® II + 7 5 <g> 7 5 



(54) 



We would like to point out here certain observation. If one considers the time ordered 
product of two \l/'s and two \l/'s then, for some choice of time arguments, one can introduce 
a complete set of (out) states obtaining for instance 

< 0\^(xi)^(x 3 )\n >< n|^(x 2 )^(x 4 )|0 > . 



The set \n > is here a Fock set of massive Schwinger bosons since these are the only 
asymptotic states. The 1-boson contribution, which leads to the appearance of a pole at 
k 2 = e 2 /n is expressible, thanks to the formulae of the LSZ formalism, through the vertex 
function. This in turn means, by virtue of the equation (pEPP , that each of the amplitudes 
above have the contributions from both k = ±1 sectors. In consequence, the polar part of G 
should, beside sectors +1 + 1 and —1 — 1, contain also traces of sectors +1 — 1 and —1 + 1. 
While calculated explicitly this leads to the substitution 7 s ® 7 s for 11 ® 1 in the last term 
of (|54"|). Why such contributions appear neither in G^ nor in remains for us presently 
unclear. 



IV. CONCLUSIONS 

The expressions for the Green functions obtained in sections [TT] and |TJ show that the 
Schwinger Model may be explicitly solved also with all topological effects taken into account. 



As argued in QlQ | , all higher functions (with no more than 4 external fermion 'legs'), thanks 
to the two Ward identities, which are preserved in the nonzero instanton sectors too, are 
the combinations of those found above. 

From the formulae ( p9|) and (^) together with (|4"2] , |5"0[ , |5"2D we can obtain, by taking the 
appropriate limits, the values of the condensates and verify the cluster property. We do not 
treat, in this work, the restoration of the cluster property as a support for the use of the 
vacuum instead of a topological one since this fact is well known and needs no verification. 
We rather show that the obtained formulae are selfconsistent and agree (in this special limit) 
with former results of other authors. 

The values of condensates 
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V = < *(x)V(x) > vac , V A =< y{x)^{x) > vac , (55) 



may be easily obtained from eq. (|39|) . If we recall that /3(0) = we get 

f %€■ 

V = -^e^ cos9 , V A = -^e- sin6 , (56) 

which agrees with the results of other authors |fl0|j28[] and constitutes the 1-instanton con- 
tribution. To verify the clusterization property let us now consider the objects defined by 



W{x,y) =< V(x)y{x)y{y)y{y) > mc , W A (x,y) =< ^(x)^(x)^(y)^(y) > r 



(57) 



The k — ±1 does not contribute here since the trace of the product of an odd number of 
gamma matrices is zero and that is what we obtain while exploiting equation flSOl). From 
the sector k = we obtain 

W (0) U y) = ^_ 2 lE+ 2K ^-e^-y)*M ? (5g) 

87r d 

•,2 



W A °\x, y) = -JL-Ji^KoW-^-vYM , (59) 
and for k = ±2 the result is 

W^- 2 \x,y) = W A +2 ''- 2 \x,y) = ^iz-^-^-yYM cos26 . (60) 



57T 3 



In particular in the limit (x — y) 2 — > — oo, where K — > 0, we have 

2 2 

W= -^e 27B (l + cos2#) , W A = -^-e 27B (-l + cos2#) , (61) 

and thanks to the identity cos 29 = 2 cos 2 $ — 1 = 1 — 2 sin 2 6* the clusterization reappears 

]AJ = V • V , = • • (62) 

We would like to devote now a few words to the 9 dependence of the Green functions we 
obtained. It is well known |7]|| that in the case of massless fermions, thanks to the chiral 
invariance of the Lagrangian, the parameter 9 may be gauged away. This happens because 
of the close connection between the chiral anomaly and the winding number. That it can 
actually be done is confirmed by our formulae (^,^,^,^). For each product of \I/ and \1/ 
in a Green function there appears a factor e _l6,7 ° in the final formulae which might obviously 
be cancelled by the appropriate chiral gauge transformation performed on the fermion field. 
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